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CN . We develop an algorithm to construct algebraic invariants for hypermatrices. We 

I then construct hyperdeterminants and exhibit a generalization of the Cay ley-Hamilton 

' theorem for hypermatrices. 

■ 1 Introduction 



- 

X 



Hypermatrices appear in several contexts in mathematics [El, |25|, E6[ and in applications such 



as in the quantum mechanics of entangled states 0, ^ , and image processing 0, |^ . Important 



;_j ■ mathematical problems in the above applications are the construction of algebraic invariants 
^ associated to hypermatrices and the determination of the minimal number of algebraically 

independent invariants. In this work we address these last two problems. 

For ordinary matrices the algebraic invariants associated to a matrix a can be obtained 
as traces of powers of the given matrix. According to the Cayley-Hamilton theorem only a 
finite number of these powers is linearly independent and therefore only a finite number of 
algebraic invariants is linearly independent. A more convenient set of invariants is given by 
the discriminants which are suitable combinations of traces and are constructed in terms of 
alternating products with the unit matrix I. 

Therefore, algebraic invariants (discriminants) can be constructed by considering all pos- 
sible products among a's and I's and these products are in a one-to-one correspondence 
with semi-magic squares of rank 2 (a semi-magic square is a square array of numbers such 
that the sum of element in each row and column gives the same result). The discriminants 
can be obtained in terms of semi-magic squares by counting all possible permutations of 



1 



indices and for practical purposes the discriminants can be constructed using a graphical 
algorithm in terms of grids which we develop here. 

We then apply the above algorithm to hyper-matrices. We restrict our considerations to 
the fourth-rank case. We obtain the corresponding discriminants, the determinant and the 
Cayley-Hamilton theorem. 

2 Matrix Calculus 
2.1 Basic Definitions 

A matrix a is a d x d square array of numbers a^j, with i,j — 1, - ■ ■ ,d. Let a and b be two 
matrices with components aij and bij, respectively. The matrix addition "+" is defined by 
the resulting matrix c = a + b with components 

Cij ciij -\- bij . (1) 
The matrix multiplication "•" is defined by the resulting matrix c = a • b with components 

d 

Cij = ^ ttik bkj ■ (2) 
fe=i 

This operation is the usual Cartesian product among rows and columns. Therefore, we have 
a ring. 

The unit element I for the matrix multiplication, a • I = I • a = a, has components Iij 
given by 

/■• = (^' fo^^ = -^' . (3) 
'■^ 1 , otherwise 

The inverse matrix a~^ is a matrix satisfying 

a-^ • a = a • a-^ = I . (4) 

In terms of componentes we have 

d d 

^ {a'^)ik akj = ttik {a~'^)kj = hj ■ (5) 
fe=i fe=i 

The conditions for the existence of the inverse a~^ will be given later on. If the inverse a~^ 
exists then a is said to be a regular matrix. 

The product of a matrix a with itself, that is a^, is the matrix with components 

d 

(a^)ii = ■ (6) 

k=l 

Products of a of an order s, a*, have the components 
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d d 

fcl = l fcs_l = l 



(7) 



s times 



s—1 times 

By definition a° = I and = a. Furthermore 



3P ■ 3?^ aP+« . 



For a matrix a the trace is given by 

d 

(a) = trace (a) = ^ 



The trace of a* is given by 



(a^)=trace(a^) = 5;](a%. 



(8) 



(9) 



(10) 



Furthermore, trace(a'') = d. 

Let u be a regular matrix. Then we can define the similarity, or affine, transformation 



a' = u ^ • a • u . 



(11) 



We can then verify that 



(a'^) = (a^) . (12) 

Therefore, the traces are invariant under similarity transformations and this set of invariants 
can be used to characterise the matrix a. 

Discriminants are a set of invariants which are given in term of traces by the relation 
among the coefficients in a Fourier expansion and the coefficients in a Taylor expansion. The 
first discriminants are given by 
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(13) 



etc. Only the first d discriminants are non-trivial, while the discriminants of an order higher 
than d are identically zero, Cs(a) = 0, for s > d. Therefore, the number of invariants 
characterising a d x d matrix a is finite and given by d. This means, by the way, that only 
the first d traces of a matrix a are linearly independent. However, discriminants are a more 
convenient set of invariants than the traces since only the first d discriminants are non trivial. 
Let us consider the derivatives dcs{si)/da.. For the first values of s we obtain 



dci (a) 

da 
5c2(a) 

da 

dcsja) 

da 

dc4{a) 



I, 

(a) I — a , 



2 
1 



(a)^ - (a^) I - (a) a + a^ 



da 3! 
etc. We can then verify that 

dcJa 



(a) -3 (a) (a^) + 2 (a' 



s-l 



{af-(a^ 



da 

Let us define the polynomials 



k=0 



a + (a) a^ - a^ , (14) 



(15) 



P,^(A) = Q(a-AI), (16) 
where A is a scalar. By means of a Taylor expansion around a, and using ([TT|), we obtain 

d 



^a(A) = E (-1)' cfc(a) A^-'^ = PfW A + (-l)'^ Q(a) . 

k=0 



(17) 



For the first values of d we obtain 



Pi(A) = A-ci(a), 

PlW = A2-ci(a)A + C2(a), 

P3(A) = A=^-ci(a)A2 + c2(a)A-C3(a), 

^i(A) = A^-ci(a)A3 + c2(a)A2-c3(a)A + C4(a). (18) 

P^(A) is known as the characteristic polynomial associated to the matrix a. We can write 
the polynomials (|I^) with a as its argument, namely P^{a) = Pi{a). Then, we obtain 

Pii^) = c,{a) a"-' = Pi-\a) ■ a + (-1)^ c,(a) I . (19) 

k=0 
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Comparison with (|T2|) shows that 



P^(a) = (-1)'^^^. (20) 



Therefore, eq. ([T6| ) can be rewritten as 



da 



a - Cd(a) I 



(21) 



The powers of a matrix a are not all linearly independent. They are related by the 
Cayley-Hamilton theorem: 

Theorem. (Cayley-Hamilton) A dxd matrix a is a root of its characteristic polynomial, 
that zs P^(a) = 0. 

Therefore, the Cayley-Hamilton theorem states that there exists a relation of the form 

• a-c,(a)I = 0. (22) 
If Cd{a) 7^ 0, then there exist an inverse matrix a^^ which is given by 

Cd(a) da 

In terms of components this inverse matrix is given by 

(a-'), = 4-%^^. (24) 

Crf(a) daij 

We have now two algorithms to compute the inverse of a matrix a. The first one is based on 
eq. (^. We have unknows {a~^)ij and equations. The condition for this equation to 
have a solution is a = det(a) 7^ 0. Then the solution is (p^). The second algorithm is based 



on the direct construction of Cd{a). If Cd{a) 7^ then the inverse matrix is given by (^31) . For 
second-rank matrices the two algorithms are equivalent and give the same result. However, 
for higher-rank matrices only the second algorithm admits a generalization. 
An additional set of relevant matrices is given by 

c,(a) = ^. (25) 

This expression is an abbreviation for the following operation. Let us consider a second 
matrix b with components bij and define c = b ■ a. Then, eq. ( ^4|) is an abbreviation for 



dh 

In terms of components we have 



b=I 
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For the first values of s the result is 



Ci a 



C2(a) = (a) a — a? 



C4(a) 
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(a)" - (a") a - (a) a' + a^ 



(a)^ - 3 (a) (a^) + 2 (a' 



1 r 



a" + (a) (2. 



We can verify that 



For s = d we obtain 



c,(a) = -(-irPr^(a) • a. 



c,(a) = -(-l)^P^-i(a) ■ a = Crf(a)I, 



where we have used (|T6|) and (jl^). Therefore 

dCdiSL 



dl 



-Q(a)I = 0, 



(29) 



(30) 



(31) 



is a statement equivalent to the Cayley-Hamilton theorem (|20|) . Let us observe that Cd(a) 
is symmetric under the interchange of a and I. Therefore, in the same way in which the 
inverse of the matrix a is defined by (|2T|) , the inverse of the matrix I is defined by 



1-^ = 1 



1 dcd{a) 



Cd(a) dl 

which is equivalent to (p9|). For the first values of d the result is 



(32) 



a — ci(a)I = 0, 

a^ — Ci(a) a + C2(a) I = 0, 

a^ — ci(a) a^ + C2(a) a — C3(a) I = 0, 

a^ — Ci(a) a^ + C2(a) a^ — C3(a) a + C4(a) I = 0. (33) 

A further set of relevant quantities is obtained as follows. Let us define 

t,(a) = - (a^) . (34) 
s 

Using ( p^ and (0) we now obtain 

. (35) 

This is the definition of the matrix multiplication in terms of discriminants which will use 
for the generalization to fourth-rank matrices. For this purpose, however, it is necessary to 
write (|33) in terms of discriminants. For the first values of s we obtain 
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ti(a) = ci(a), 

t2(a) = -c^(a) - C2(a) , 

t3(a) = ic5(a) -ci(a)c2(a) + C3(a), 

t4(a) = ici(a) -Ci(a)c2(a) + ci(a)c3(a) + ^C2(a) -C4(a). (36) 
2.2 Index Notation 

For the purposes of generalisations to higher-rank matrices the notation above is not ad- 
equate. We have found that an index notation similar to that of tensor analysis is more 
convenient. In this case is necessary to distinguish between covariant and contravariant in- 
dices. According to this scheme, the matrix a of the section above will be now a second-rank 
covariant matrix aij. Analogously, a second-rank contravariant matrix b has the components 

In order to obtain the index version of the matrix multiplication for covariant matrices, 
let us consider a second-rank contravariant matrix I with components 

/*■? = /"'"' if ^ ~ i' /^Y) 
1 , otherwise. 

Then the matrix multiplication is given by 

dj = ttik I^^ kj . (38) 

(The summation convention over repeated indices is assumed.) In this sense, I behaves as 
a metric. With these conventions we can reproduce all the definitions and results of the 
previous section, but now product and traces are defined and constructed with P^. 

The similarity transformations are now constructed in terms of a matrix u of mixed 
covariance with components u^j and the inverse with components {u'^Yj, such that 



u^,iu-'r^ = {u-J,u', = 5]. (39) 



Then 



a' = u ■ a • u 



b' = u"^ ■ b ■ , (40) 



or, in terms of components. 



{a')ij = ttki u^i u'-j , 

{br = b'^{u-')\{u-y,. (41) 

Then, the discriminants (|T2p are invariant under the similarity transformations 



7 



2.3 Alternating Products and Discriminants 

In terms of the components of the unit matrix I we can define the following symbols 

ghjl - isja _ jiijl . . . jisja _ J_ . . . jisjsjl ^^2) 

where Pj denotes the sum over all permutations with respect to the second index j; |[- ■ ■]! 
denotes complete antisymmetry with respect to the indices j's or, equivalently, with respect 
to the indices i's. Due to the antisymmetry the symbols q are non trivial only for s < d. 
For the first values of s q is given by 

^«lil*2j2i3j3 _ J_ I'jilil p2j2 ji-ij-i _ ^jnjl p2jz p3j2 _j_ jhjs p2j2 psjl 

_|_j«u'2 p2ji jiaja^ _|_ ^pij2 p23z jhji _j_ jnis p23i j«3i2^j (43) 
Then, the discriminants are given by 

Cs(a) = q^^'-'^^' ai,,, ■ ■ ■ a^^j^ . (44) 



Let us observe that 



Therefore 



jh---idji---jd _ j|[nji ... pdjd]\ — ^h---id ^h-id ^45^ 



Q(a) = -e^--^'*e^--^''a,,,, (46) 



In matrix calculus the usual definition of the determinant is given by 



det(a) = le^--''6^--^'^a,,,, ■■■a,,,,. (47) 



Therefore 



Crf(a) = det(a) . (4J 
Let a = det(a). If a 7^ 0, then eq. (^) can be rewritten as 



In terms of components 



Therefore 



a- = i^. (49) 

a osl 



«« = i|^. (50) 
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7«J = ]. 1 Jh-i{d~i) Jh-i{d-i) r, „. 

{d-\)\a '^Hd-i):'(d-i) 



(51) 



We can verify that 



(52) 



The relations (pUHsTI) are the usual definitions of the inverse of a matrix in matrix calculus. 



2.4 Semi— Magic Squares 

Monomial algebraic invariants are characterised by the way in which the indices of a are 
contracted with the indices of I. This contraction scheme can be represented by a semi- 
magic square. In order to construct algebraic invariants we consider products of a's and I's. 
In order for the result to be an invariant all indices must be contracted. This means that we 
must consider an equal number n of a's and I's. Since both a and I have two indices each, a 
can be contracted at most with 2 indices belonging to I's and the same is true for I's. The 
way in which the n a's are contracted with the n I's can be represented by an n x n square 
array of numbers s, where the components s/j denote the number of contractions between 
the Jth a and the Jth I. Graphically 

Il ■ In 

ai / sii ■ ■ ■ si„ \ 

s = ; ; ; . (53) 

aji y "Snl ■ ■ ■ Snn j 

If the /th a is contracted once with the Jth I, then we write sjj = 1; If the Ith a is contracted 
twice with the Jth I, then we write s/j = 2; If there is no contraction between the Ith a 
and the Jth I, then we write s/j = 0. Since all indices must be contracted the sum of the 
elements of each row and each column must be equal to 2, that is, 



7=1 



= 2. (54) 

j=i 



Arrays with this property are known as semi-magic squares |]rB|. Therefore, the number of 
possible algebraic invariants is determined by the number iJ„(2) of semi-magic squares s. 
Semi-magic squares of rank r are defined by the relations 



1=1 

n 



j=i 
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The number of different possible semi-magic squares Hn{r) is given by [0, |T^, |T8 



HAr) 
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(56) 



wfiere i? = (r + l)(r + 2)/2. 

For r = 2 the resuh is: Hi{2) = 1, 1/2(2) = 3, H-i{2) = 21, 1/4(2) = 282. For n = 1 and 
n = 2 the corresponding semi-magic squares are given by 
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(57) 



etc. Since each column and each row represent the same matrix, semi-magic squares which 
are related by the permutation of rows and/or columns represent the same algebraic in- 
variant. Therefore, semi-magic squares can be classifed into equivalence classes related by 
permutations of rows and/or columns. Therefore, we need to take care only of the repre- 
sentatives Sr^n = = 1, ■ ■ ■ ,Pr{n)} foT cach equivalence class, where Pr{n) is the number 
of equivalence classes for rank r and order n. For r = 2 the number of equivalence classes 
P2{'n) is given by the number of integer partitions of n, taht is, p{n). For the first values of 
n the representatives of each equivalence class are 
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VO 1 1/ 

The algebraic invariants which each semi-magic square represents are given by 
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Let us observe that block semi-magic squares can be decomposed in terms of lower order 
semi-magic squares as 
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2.5 Permutations 



Each algebraic invariant is represented uniquely by a semi-magic square. Therefore, semi- 
magic squares are adequate to represent the different algebraic invariants. A discriminant of 
order n is a linear combinations of the semi-magic squares of order n. In order to determine 
the coefficients of this linear combination we proceed as follows. 

A matrix a can be contracted with a second matrix a in a number of ways which depend 
on the possible permutations, 7r(n), of the indices in the second matrix. For each possible 
permutation we obtain a signed semi-meagic square /cj, i = 1, • • • , 7r(n). The discriminant is 
then given by the general formula 



Cn(a) 



i=\ 



(61) 



where 7r(n) is the number of possible outcomes ki after considering all possible permutations. 
The numerical factor 1/n! in front of the previous expression is conventional. 

The signed semi-magic squares hi are obtained as follows. For n — 1 the possible per- 
mutations are 



TT 




(62) 



where the subindex denotes the parity of the given permutation; Therefore, 7r(2) = 2. We 
can always choose to keep fixed the indices of the first matrix and represent it by a square 
with I's as entries in the diagonal. For n = 2 we have 



TTo = 



(63) 



The possible outcomes are 
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Then, the discriminant is given by 
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For n — 3 there are 6 possible permutations, 7r(3) = 6, given by 
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The discriminant is given by 
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For n = 4 there are 24 possible permutations, 7r(4) = 24, namely. 
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The discriminant is given by 
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The Caylcy-Hamilton theorem is obtained by deriving the corresponding discriminant 
with respect to I. In the language of semi-magic squares this means that one of the I's, 
one at a time, is not present. Therefore, the indices of that column are not contracted. We 
indicate that derivative by the same original semi-magic square and 'resalt' the uncontracted 
column by boldface numbers. The resulting Cayley-Hamilton theorem is now written as 
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where [2] stands for the unit matrix I^-*^. Contracting with I, [2] = d, we obtain (|5^, ( p3| ) 
and 



+ 



1/ 



C4(a)[2] = 0, (71) 



2.6 Graphical Construction of Invariants 

Each algebraic invariant is represented uniquely by a semi-magic square. On the other 
hand, semi-magic squares are obtained by considering all possible permutations of indices. 
However, for large values of n this algorithm becomes unpractical. In order to avoid this 
difficulty we now develop a graphical algorithm for the construction and characterisation of 
algebraic invariants which allows to simplify this task. Let us represent the matrix a by a 
vertical grid with two boxes, namely 



- (72) 
The product of n matrices is represented by 



... (73) 

12 n 

Each algebraic invariant is characterised by the way in which indices are contracted. We 
can always choose to keep fix the indices of the first row and look at how the indices in the 
second row are contracted with the indices in the first row. A void grid indicates that no 
permutation has been performed 



(74) 
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A permutation of the indices ith and jth is indicated by 



A cyclic permutation is indicated by 



(75) 
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A double permutation is indicated by 
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In this case, however, it is necessary to take into account the sense in which the permutation 
is performed. There are 2 possibilities 



(7J 



When the sense of the permutation is irrelevant we use the right-oriented grid. The next 
possibility is 



(79) 



In this case there are 6 possibles senses for the permutation, namely. 



o- 


























(» 
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-o-> 
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'1 > 



















CM 




-o 





















(80) 



As in the previous case, when the sense of the permutation is irrelevant we use only the 
right-oriented grid. 

The semi-magic square corresponding to a given grid is obtained as follows. The number 
of empty boxes in each column corresponds to the diagonal entries in the semi-magic square. 
The lines correspond to the off-diagonal terms. For example 

/2 0\ 

Oil. (81) 

The multiplicity is the number of possible ways in which the given permutation can be 
performed over the n boxes (indices). The parity is given by the number of lines for the 
permutation. For example 
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(82) 



Let us perform the explicit construction of the semi-magic squares and discriminants for 
the first values of n. For n = 2 the result is 



2 
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1 1 
1 1 



(83) 



Therefore, we obtain (1621). For n = 3 the result is 
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Therefore, we obtain 



. For n = A the result is 
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(85) 



-6 



/I 1 ON^ 

110 

11 

VI 1/ 



Therefore, we obtain (|66D. 

The number of semi-magic squares depends on the number of possible permutations. For 
large values of n and r this counting becomes quite involved. Therefore, it is advisable to 
have an easy recipe to obtain the correct counting. 

For each order we have a different number of possible permutations -P„. We can represent 
them as 
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P2 = 0-1, 

P3 = 0-3-1 + 2-2, 

P4 = O-e-l + S-l^ + 8- 2- 6- 3, 

P5 = 0-10- 1 + 15- 1^ + 20 -2 - 20- (21) -30- 3 + 24 -4, 
Pg = - 15 - 1 + 45 - 1^ - 15 - 1^ + 40 - 2 - 120 - (2 1) 

+40 - 2^ - 90 - 3 + 90 - (3 1) + 144 - 4 - 120 - 5 . (86) 

The number of terms involving some given permutations is given by the coefficients in (|5BD. 
This can be easily verified in the graphical construction above. 



3 Hyper mat rices. The Fourth— Rank Case 

There is not a natural multiplication operation for hypermatrices in the sense that the 
product of two hypermatrices be again a hypermatrix of the same rank. Therefore, the 
construction of algebraic invariants must be performed in a different way which, by the way, 
gives a clue for the definition of amultiplication operation for hypermatrices. We now turn 
to the construction of algebraic invariants for fourth-rank matrices. This algorithm can be 
easily extended to matrices of an arbitrary even-rank r; for the odd-rank case see section 4. 

Our construction is based on alternating products. To this purpose let us consider a 
fourth-rank matrix A with components Aijki and a fourth-rank matrix A with components 
A*-''^' as a fourth-rank generalization of the unit matrix. For practical purposes the discrim- 
inants are better represented by semi-magic squares of rank 4. They are constructed using 
the graphical algorithm of the section above. 



3.1 Alternating Products and Discriminants 

Let us now introduce a fourth-rank unit matrix A*-'^'. In order to fix the ideas we can 
consider a unit matrix defined by the relations 

^^Jkl^n, iU=J = k = l, 

I , otherwise. 



However, other definitions are possible and they are currently under study |2^. In a way 
similar to (Eol) we define 



Qiljlk\h---iajaksla _ }_ . . . ^isjaksls]\ 

" s\ 
For the first values of s we obtain 



Qi 



ijkl 



Q 



i\j\k\hi2j2k2l2 



^ijkl 
1 



^^hjikih ^i2j2k2l2 
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_^^iijikil2 ^i2j2k2h _j_ ^iijik2li ^12^2^1/2 
_^^iij2kih ^i2jik2l2 _|_ ^i2jikili ^iij2k2l2'j 

_^^^iijik2l2 ^i2j2kili _|_ ^iij2kil2 ^i2jlk2h _j_ ^«2iifei'2 ^ni2fc2'l^j (89) 

etc. Then, the discriminants are defined by 

Cs{A) = Ql^nkrh-.sM A^^^^^^i^ . . . ■ (90) 

Let us now consider a mixed second-rank matrix U with components Wj and its inverse 
U^^ with components {U'^Yj. Then we define a similarity transformation for fourth-rank 
covariant and contravariant matrices as 



A' = [A, U, U,U, U], 

B' = [B,V-\V-\V~\V-'], (91) 



defined by the components 



( A'\. . . . — A. . . . rrji. rrj^. rT^'>. tf^*- 

^^/yii2i:ii4 ^ B^^^^^''^'' {u~y^j^{U'y^j^{U'^y''j,^{U'^y^j^. (92) 

Then, the discriminants are invariants under this kind of similarity transformations. 

The determinant for a higher-rank matrix can be defined in complete analogy with the 
definition for ordinary matrices. Therefore, in analogy with (^) we define 

A = det(A) = Cd{A) . (93) 
Let us observe that a relation similar to (^) holds for fourth-rank symbols, namely. 



Therefore 



Qh-id3i-3dk\-kdh--la _ J_ ^| [njifci/i ^idjdkdldW _ }_^^h-id ... ^h-ld ^g^^ 



A = det(A) = ^ e'--'^ ■ ■ ■ e'-'''* ■ ■ ■ Am ■ (95) 



In analogy with (47) we define 



A dAijki 
Then 



A^jki = 1 . (96) 



1 1 

{d-iy. A 

This matrix satisfies 



A''"' = T^-^ - ■ ■ ■ ^'''-'''"^ Annk.1. ■ ■ ■ A,,_,,3,d-.,k,,_,,kd-r^ ' (9^) 



A'^-^^^' Ak,k2k, = 5] . (98) 
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The definitions ( p2[ ) and (p3| ) were used in previous works [0, |2y, |21], ^ concerning the 
apphcations of fourth-rank geometry to the formulation of an alternative theory for the 
gravitational field. 

Let us observe that Cd{A) is symmetric under the interchange of A and A. Therefore, 
in analogy with ( pOj ) the inverse of the matrix A is defined by 

a,(A] OA ■ 



In terms of components 



This matrix satisfies 



In terms of components 



_ 1 dC,{A) 
~ CjA) ■ 

A-i ■ A = A ■ A-i = I. (101) 
Aik,k,k,A^'''''' = Si. (102) 



The eq. (^J) can be rewritten as 



CdiA)A = 0. (103) 



This is the statement of the Cayley-Hamilton theorem for hypermatrices. 

As an example of the relation above let us consider the simple case d = 2. The determi- 
nant (|92|) is then given by 

A = ^1111 ^2222 — (^1112 ^2221 + ^1121 ^2212 + ^1211 ^2122 + ^2111 ^1222) 

+ (^1122 ^2211 + ^1212 ^2121 + ^2112 ^122l) ■ (104) 

The components of the matrix A*-''^' are given by 

/lllll _ 1 A /11112 _ 1 A /11121 _ 1 A /11211 _ 1 A 

^ —^^2222, ^ — — ^4^2221, ^ — " ^4 ^2212 , ^ ——^^2122, 

/12111 _ 1/1 /III22 _ 1 /I /|1212 _ 1 /I \1\V2 _ 1 /I (105) 

^ ——^^1222, ^ —^^2211, ^ —^^2121, ^ —^^1221, 



and similar expressions for the other components. In order to check the validity of eq. (|9 
let us consider the cases (11) and (12). We can then verify taht 

A^'^'^A^,,^ = 0, (106) 
and similar relations for the other indices. 
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3.2 Semi— Magic Squares 



The algebraic invariants which can be constructed in this case are given by the semi-magic 

squares of rank 4. Their number is Hi{A) = 1, i?2(4) = 5, H^{A) = 120, i?4(4) = 7558. As 
for the second-rank case we must take care only of the representatives for each equivalence 
class. The number of equivalence classes P4(n) is given by the generating function 



1 



n=0 



For the first values of n P4{n) is given by 



P4{n) = {1, 1, 3, 9, 36, 
The representatives for each equivalence class are 
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Let us observe that in (|106|) there are 40 magic squares of order 4 rather than the 36 imphed 
by (|TU§; the couples (19,20), (25,26), (34,35) and (36,37), labelled by a, b, c and c, 
respectively, have the same connectivity structure but, since rows and columns corresponds 
to different matrices, they represent different invariants. 



3.3 Construction of Discriminants 

Each semi-magic square represents one and only one algebraic invariant. However, in this 
case, the corresponding invariant can no longer be represented by mean of traces, as was 



done for ordinary matrices. The semi-magic squares of order 1 and 2 in ( 104 ) correspond to 
the following invariants 



(4) 



ijkl 1 



4 
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Ahjlkih A, . . , \i2j2k2l2 /I, . . , 
^ ^h^2J2k2 ^ ^t2nji«:i ' 
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2 2 
2 2 



k2h 



^k2l2iiji 



(110) 



It is obvious from the expressions above that semi-magic squares are more practical for 
representing algebraic invariants. 

The corresponding discriminants are linear combinations of the monomial algebraic in- 
variants (semi-magic squares) of the same order. In order to determine the coefficients of 
this linear combination wc proceed in a way similar to that for ordinary matrices. The 
matrices can be contracted according to the allowed number of possible permutations. The 
possible permutations are the same as described in sec. 3.2. However, this time they must 
be summed three times. It is obvious that the number of terms which must be computed 
growth very fast, as (n!)^~^. Therefore, even when this algorithm provides a direct answer, 
a more practical way to evaluate the coefficients is necessary. Let us therefore consider the 
graphical algorithm developed in sec. 2.6. 

For n = 1 there is only one possibility, namely 



Ci(A) = (4) 



For n = 2 we have 
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The coefficients in this expressions are obtained from 

P| = (0-l)'^ = 000-3001 + 3011-lll. 
The resulting discriminant is 

C2{A) 
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(113) 



(114) 



For n — 3 the grids are given in Appendix A. We can verify that the coefficientes are 
given by 



P| = (0 - 3 ■ 1 + 2 ■ 2)^ 

= 000-9001 + 27 011 + 600-2- 27 111 

-36-0-1-2 + 54-1-1-2 + 12-0-2-2-36-1-2-2 + 8- 2- 2- 2. (115) 

The discriminant is 
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For n — 4: the grids are given in Appendix B. Now we do not verify that the coefRcientes 
in the relation 



p3 ^ (0-6- 1 + 3- 1^ + 8- 2-6-3)^ 
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corresponds to that given in the Appendix B. Instead, we use them to check the correctness 
of the coefficients there. The discriminant is 
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3.4 The Cayley— Hamilton Theorem 

Let us write the Cayley-Hamilton theorem in terms of almost-magic rectangles. For the 
first values of d we obtain 
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where [4] stands for the unit matrix A^^. Contracting with A, [4] = d, we obtain ( |113| ), ( p.l4| ) 
and ( |115| ), respectively. Therefore, the above is the Cayley-Hamilton theorem for fourth- 
rank matrices. The expression of the Cayley-Hamilton theorem for d = 4 is exhibited in the 
Appendix C. 

In order to ilustrate our result let us consider the simple case of a completely symmetric 
fourth-rank matrix Gijki in two dimensions, j, k,l = 1, 2. Then, the determinant, according 



to (110) and ( 113 ), is given by 



C2(G) — Gllll 6*2222 ~ 4 6*1112 6*1222 + 3 G*^;^22 



(120) 



In terms of G's and A's the almost magic rectangles in ( |1 16| ) are given by 



4 

4 , 

^ abed 

3 l^ 

1 3 / abed 

2 2\ 

2 2 / ^f^^ 



(A*-'^' Gijkl) Gabcd 5 

GA ijkl 
{a\ijk ^ ^l\bcd) , 

GA ijkl /~i 
{ab\ij ^ ^kl\cd) ■ 



(121) 



In detail we have 



4 
4 



1111 



1111 + 6*2222) G 



1111 ; 
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4 

4 

4 

4 

3 1 

1 "^/iiii 
3 1\ 1 



1112 



1122 



(Cliii + 6*2222) ^1112 , 

{Gnu + (^2222) 6*1122 , 

^1111 + 6*1112 G1222 ; 



^ 0/ 7 (4 Giiii 6*1112 + 6*1112 G2222 + 3 G1122 6*1222) 

'^^1112 4 

31\ If/^ . ^ .^2 I ^2 



-, q / 9 [(6*1111 + G2222) 6*1122 + + ^1222] 5 

"^/1122 ^ 



2 2 
2 2 



1111 



(<^2 I ^2 

"-^1111 "T '-^1122 ' 

6*1111 6*1112 + 6*1122 6*1222 ; 

[(6*1111 + 6*2222) 6*1122 + 2 6*^112 + 2 6*^222] ' (122) 
and similar expressions for the other indices. Then we have 



2 2 

2 2\ _ 1 
2 2/^^22 3 



")„„-' (1 Dun^'il 2)„,r^^<^>M- ^ °- 
^0 4)„,r'(? 3)„,,-K2 2)„,r^^<^>W'- ^ 

and similar results for the other indices. 

3.5 Generalised Multiplication and Ring Structure 

Now we introduce a definition for the multiplication of fourth-rank matrices which allows to 
give a ring structure to the space of fourth-rank matrices. Such as we defined determinants 
by means of an inverse analogy, we define the product of fourth-rank hypermatrices with a 
similar relation. To this purpose, in analogy with (p3|), we define 



Til 


:a) 


= C^i(A), 








T2I 


:a) 


= 


- C2{A) , 






T3I 


:a) 


= IcUA) 


-6*1 (A) 6*2 (A) 


+ 6*3( 


;a). 


T4I 


:a) 


= IctiA) 


-C2(A)6'2(A) 
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Now, we invert the relation (^5[) to define powers of fourth-rank matrices as 
For the first values of s we obtain 



(A) 



-12 



+12 
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0^ 
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(3 
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2 






2 
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1 




2 
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/3 1' 
3 1 




3 1 0" 
12 I 1 2 1 
1 3 



The expression for (A^) is displayed in Appendix D. 
Now, the discriminants can be written as 



12 



/3 1 
1 1 
\0 2 





(2 


2 






12 
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1 


I) 






[o 


1 







where 



Co(A) 

Ci{A) 

C2iA) 
CsiA) 



1, 

[A] , 
1 r 

2 
1 

3! 
1 

4! 



Al 



\A] 



A^ 



3 [Al 



A^ 



[Al^ - 6 \A? 



A^ 



A^ 



+ 8 fAl 



A^ 



+ 3 



A^ 



[A] = Trace ((A)) = A ■ (A) . 
The Cayley-Hamilton theorem can now be rewritten as 



A^ 



(A)-Ci(A)I = 0, 

(A2)-Ci(A)(A) + C2(A)I = 0, 

(A3)-Ci(A)(A2) + C2(A)(A)-C3(A)I = 0, 

(A^)-Ci(A)(A-^) + C2(A)(A2)-C3(A)(A) + C4(A)I = 0. 

which is formally exactly the same as (p3D. 
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Since I is a fourth-rank matrix, the multiphcation defined above is not a binary operation. 
In fact, we can define binary, trinary and tetranary products as 



(A, A)^ = AIA, 
(A, A, A)3g = AIAIA, 
(A, A, A, A)2i7 = AI AIAIA. (130) 

The subscript denotes the number of terms involved in each expression. Then, the products 
defined above are given by 



(A^) = (A, A), 

(A3) = (A, A, A) + (a^, a) , 

(A^) = (A, A, A, A) + (a^, a, a) + (A^ a) + (A^ A^) . (131) 
Let us finally observe that 

(A*') ■ {A") ^ (AP+'?) . (132) 

4 The Third-Rank Case 

In this section we study the third-rank representative of odd-rank hypermatrices; 

all other odd-rank cases can be dealt with in a similar form. 

Let us start by considering a third-rank matrix A with components Aijk and a third- 
rank matrix A with components A^^^ as a third-rank generalization of the unit matrix. For 
practical purposes the discriminants are better represented by semi-magic squares of rank 
3. They are constructed using the graphical algorithm of the section above. 

Let us introduce a third-rank unit matrix A*-''^. In order to fix the ideas we can consider 
a unit matrix defined by the relations 

^ijk _ 1 1, if ^ = = _ (133) 
L U , otherwise. 

However, other definitions are possible and they are currently under study PB| . 

Let us therefore proceed to the construction of discriminants for third-rank matrices. 
For n = 1 there is only one possibility, namely 



For n = 2 we have 



Ci{A) = (3) 



(134) 
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1 2 

2 1 



(135) 
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The coefficients in ( |135| ) are obtained from 



p| = (0 - 1)2 = ■ - 2 • • 1 + 1 ■ 1 

The corresponding discriminant is given by 



3 
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2 1 
1 2 



For n = 3 the corresponding grids are given by 
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(136) 



(137) 



(138) 



The coefficients are obtained from 



= (0-3- 1 + 2 -2) 
= 00-601 

Then, the discriminant is given by 



9-1-1 + 4-0-2-12-1-2 + 4-2-2 



(139) 




(140) 



For n = 4 the grids are given in the Appendix C. The coefficients are obtained from 
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(0-6- 1 + 3- 1^ + 8- 2-6 -3)2 
O O+ O - l+ l - l+ O - 12+ 0-2 



+ ■ 1- r +1-2+ 0-3 4 
+ ■ 1^ • 2 + -2 • 2 + -l^ ■ 3 



■1 • 3 + 1^ ■ V 
■2 -3 + -3 -3. 



The discriminant is then given by 



1 

24 



-24 



+3 



+36 



^3 ON^ 
3 
3 
VO 3/ 
^3 ON^ 
2 10 
111 
VO 1 2/ 
/2 1 ON^ 
12 
2 1 
VO 1 2/ 
^2 1 ON^ 
2 1 
10 11 
VO 1 1 1/ 



+ 



/3 ON^ 
3 
2 1 
VO 1 2/ 
^3 ON^ 
111 
111 
VO 1 1 1/ 
/2 ^ 
2 1 
10 2 
VO 1 1 1/ 
/2 1 ON^ 
10 11 
111 
VO 1 1 1/ 



+ 24 



-6 



24 



24 



^3 ON^ 
2 10 
2 1 
VO 1 2/ 
^2 1 ON^ 
2 10 
2 1 
VI 2/ 
/2 IN^ 
2 1 
2 1 
VI 1 1 0/ 
^1110^ 
111 
10 11 

VI 1 1/ 



+ 12 



(141) 



• (142) 



Let us now consider the previous expressions for d = 2 for a third-rank hypermatrix a 
with components aijk- Then, the relevant discriminants are 



— C^OOO + (^111 5 
C2(a) — C^OOO C^lll ~ - (c^OOl Alio + C^OIO filOl + c^ioo floii) • 

According with (^) the components of the inverse hypermatrix are given by 

1 9C2(a) 



ijk 



In detail we obtain 



Co (a) da 



ijk 



(143) 



(144) 



,000 



■ am 



,001 



3C2(a) 



^110 , 
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3C2(a) — ' 



and similar relations for the other components. Then, we obtain 

ao,ja^'^ = 1, 

ao^,a'*^■ ^ 0, (146) 
and similar relations for the other components. Therefore 

aMa^''^6{. (147) 

Therefore, even when C2(a) is an invariant for a third-rank hypermatrix it does not have 
one of the essential properties of a determinant. 



In order to understand this result let us consider the relations similar to (^0]) . We define 



s\ 

For the first values of s we obtain 



Qiljiki-isjsks _ ... ^isjaksW (148) 



_(^/\hjik-L I\i232k2 _|_ I\i\hk2 ^i232ki 
_|_^ni2fei ^i2jlk2 _j_ ^i2jiki ^ni2fc2^ 

_|_^^nilfe2 ^«2j2fcl _j_ ^«lJ2fcl ^«2jlfc2 _j_ ^«2jlfcl ^ni2fc2^j (149) 

etc. Then, the discriminants are defined by 

Cs{A) = Ql^nkr-.s3sk. A^^^^^^ . . . A^^^^^^ (150) 

Let us now consider a mixed second-rank matrix U with components U^j and its inverse 
U^^ with components {U'^Yj. Then we define a similarity transformation for fourth-rank 
covariant and contravariant matrices as 



A' = [A, U, U, U, U] 



B' 



= [B, U-\U-i,U"i, U-i], (151) 



defined by the components 



J3 . 
13 ; 



( A'\. . . — A. . . rrn. ttJ2, rrja 

\^ Al«2«3 ^J1J2J3 ^ «1 ^ «2 ^ 

^Qyl^2i3 ^ ^hhh {U~^Y^.^ (^"^)''j2 (^"^)'%3 • (152) 
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Then, the discriminants are invariants under this kind of similarity transformations. 



Let us observe that a relation similar to (|43|) does not holds for third-rank symbols, 
namely, 

1 1 

Qii-idji-jdki - kd _ _^|[njiA:i . . . ^idOdkdW _^ _ ^'^i-'^d . . . ^ki- ka (153) 

The determinant for a higher-rank matrix could be defined in complete analogy with the 
definition for ordinary matrices. Therefore, Crf(A) 7^ det(A). In fact 



C,{A) = Qy^^n-,dk.-k, ^^^^^^^ . . . ^^^^^^^ ^ ^ . . . ^^^^^^^ . . . ^ 



idjdkd ■ 

(154) 
However 

A = det(A) = 1 e^--^" ■ ■ ■ e^--^^ A,,,,,, ■ ■ ■ A,,,,^, ^ . (155) 
This result holds for any odd-rank hypermatrix in odd dimension. In fact 

det(A) = e-"--^^ . e'r.-.rd^ A,,,...,^, ■ ■ ■ A,,,...,^, ^ . (156) 
7 ^ ' 

^t™*''^ d times 

This result is due to the odd number of e's which add all contributions to zero. This result 
can be easily verified for the simple case at hand, r = 3, = 2. We have 

g 1 2 giii2 g 1 '2 Ai-^j-^i^-^ ^i2j2k2 = • (157) 

In order to obtain some indication as to the correct way to define a determinant for odd- 
rank matrices, let us consider the simple case of a completely symmetric third-rank matrix 
a with components ttijk in dimension d. Then, let us look for a solution a*-'^ to the equation 

a''^'"' a,k,k2 = S] , (158) 

which defines an inverse matrix in a way similar to (0). The number of unknowns is 
d{d + l){d + 2)/6 while the number of equations ( |ll(j| ) is d"^. This algebraic system of 
equations is underdetermined, except for d = 2. In this last case the solution is given by 



— (aooo c^iii + 2 aQi;^ — 3 aooi ctoii 0,111 
a ^ ' 

— I Z 0-111 Q-ooi '3'000 C^Oll (^111 C^OOl '3-011 

(159) 



and similar relations for a^^^ and a°^^, where 



'^22 3/i322 /\ 

'^000 '^111 '^000 '^001 C^Oll 0-111 + 4 

c^ooo ctoii + '^111 '^001 ~ 3 ctooi '^011 ■ (160) 



1 

Let us observe that 
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,000 



,001 



1 da 
a daQQQ 
1 1 da 



or 



ijk 



3 a (9aooi 



1 da 



a da. 



ijk 



fl611 



(162) 



Therefore, the role of the determinant, in a way similar to (^71) and (@), is played by a. 

Let us observe that for any matrix A the product of e's and A's have sense only for an 
expression which is at least quadratic in A's and A's must be of even rank. On the other 
hand, a is quartic in a's and in order to relate a with some determinant it must be with 
the determinant of some hypermatrix A which is quadratic in a's, therefore a sixth-rank 
matrix. Let us therefore consider a sixthh-rank matrix A with components Ai^ 



Since 

the rank is even we can use an extension of the definition (pO| ) to sixth-rank matrices. For 
d = 2 we obtain 



A=- e''^' 
2 



A- .4 



(163) 



Let us therefore construct discriminants for rank 6 hypermatrices. For n = 2 the corre- 
sponding grids are 
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4 2 
2 4 



1 5 
5 1 



(164) 



The coefficientes are obtained from 



P| = (0-1)^ 

= OOOOO-500001 + lOOOOll 

-lOOOlll + 501111-lllll. (165) 

Then, the discriminant is given by 
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C2{A) = - 



6 



+ 15 



10 



(166) 



For n — 3 the grids and semi-magic squares are given in Appendix D. The coefficients 
are obtained from 



p5 
-^3 



(0 - 3 • 1 + 2 • 2)^ 

OOOOO-15OOOO1 + 90OOO11 
+10 • - 2 - 270 0111- 120 01-2 
+405 •01111 + 540 0011-2 + 40 000-2-2 
-243 • 1 • 1 • 1 • 1 • 1 - 1080 •0111-2-360 001-2-2 
+810 • 1 • 1 • 1 • 1 • 2 + 1080 •011-2-2 + 80 00-2-2-2 
-1080 •1-1-1-2-2-480-0-1-2-2-2 + 720-1-1-2-2-2 
+80 •0-2-2-2-2-240-1-2-2-2-2 + 32-2-2-2-2-2. 



(167) 



The discriminant is given by 



C3(A) 



1 
6 




(168) 
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Let us now return to the two-dimensional case. Then we have 

^ = 2 [^000000 ^111111 — (^000001 ^111110 + ^000010 ^111101 + ^000100 ^111011 

+^001000 ^110111 + ^010000 ^101111 + ^100000 ^Olllll) + (^000011 ^111100 
+^000101 ^111010 + ^001001 ^110110 + ^010001 ^101110 + ^100001 ^011110 
+^000110 ^111001 + ^001010 ^110101 + ^010010 ^101101 + ^100010 ^011101 
+^001100 ^110011 + ^010100 ^101011 + ^100100 ^011011 + ^011000 ^100111 
+^101000 ^010111 + ^110000 ^OOllll) ~ (^000111 ^111000 + ^001011 ^110100 
+^010011 ^101100 + ^100011 ^011100 + ^001101 ^110010 + ^010101 ^101010 
+ ^100101 ^011010 + ^011001 ^100110 + ^101001 ^010110 + ^110001 ^ooiiio)] • (169) 

Let us consider a third-rank matrix a with components aijk- Then, we define the sixth-rank 
matrix 



4 • ■ , 



jikii2j2k2 



^iljlki ® '-^12^2^2 

T i^hjiki ^1232^2 + ^iijik2 ^i2j2ki + ^iij2k\ ^i2j\k2 + ^i2jikx (^iij2k2) ■ (I'^O) 



For d = 2, A has 64 components. We exhibit 32 of them; the other 32 components are 
obtained by interchange of O's and I's. We obtain 



^000000 


2 

~ ^000 5 








^000001 


= ^000010 = 


^000100 — 


^001000 


— ^010000 — ^100000 — flooo flool 


^000011 


= ^010001 = 


^001010 = 


^011000 


= 2 ('^000 c^oii + c^ooi c^oio) ; 


^000101 


= ^100001 = 


^001100 = 


^101000 


= 2 ('^000 c^ioi + c^ooi ctioo) , 


^000110 


= ^100010 = 


^010100 = 


^110000 


= 2 ('^000 Q-iio + c^oio ctloo) , 


^001001 


2 

~ '^001 ' 








^010010 


2 

~ '^010 ' 








^100100 


2 

~ '^100 > 








^000111 


= ^100011 = 


^010101 = 


^011100 


= ^110001 




= ^ (c^ooo c^i 


11 + ^001 0- 


110 + C^OIO C^lOl + (^100 Ooll) , 


^001011 


= ^011001 = 


C^OOl C^Oll ) 






^001101 


= ^101001 = 


O-OOl Oloi , 






^001110 


= ^110001 = 


c^ooi c^iio ■ 







(171) 



The determinant, using (|124|) , is given by 
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A = \c, (172) 



where 



_ ( 2 2,2 2,2 2,2 2\ 
— 1*^000 '^111 + '^001 '^110 + '^010 '^101 + '^100 '^011 J 

~2 [aooo C'lii (c^ooi c^iio + c^oio c^ioi + c^ioo c^oii) 

+aooi floio c^ioi o-iw + c^ooi c^oii c^iio o-im + c^oio c^oii c^ioi Oioo] 

+4 (aooo c^oii c^ioi c^iio + c^ooi <^oio <^ioo c^iii) , (l'^3) 

is the determinant introduced by Cayley |p. 

Let su mention that there is a further algorithm to construct the Cayley hyperdetermi- 
nant. Let us consider the symmetric matrix 

Qij = 0'ikik2 o-jhh ^^^^^ ^^^^^ ■ (l'^4) 

Then we have 



Qoo = 2 (aooo ^oii — ctooi aoio) , 

Qoi = aooo am — aooi ano — aoio aioi + aon aioo 

Qii = 2 (am aioo — aiio aioi) . 



(175) 



The determinant of this matrix is given by 



9oi 



det(g) = goo9n 
4 (aooo aoii — aooi aoio) (am aioo ~ ano aioi) 
— (aooo am — aooi ano ~ aoio aioi + aon aioo)^ 



-C. 



(176) 



This algorithm appeared in 



5 Concluding Remarks 

We have developed an algorithm to construct algebraic invariants for hypermatrices. We 
constructed hyperdeterminants and exhibit an extension of the Cayley-Hamilton theorem 
to hypermatrices. 

These algebraic invariants were considered by Cayley see Q for an updated account. 

Higher-rank tensors look similar to hypermatrices and the results presented here are a 
first step for the construction of algebraic invariants for higher-rank tensors. Higher-rank 
tensors appear in several contexts such as in Finsler geometry H, |T6[, fourth-rank gravity 



19, pn, 21, 22], dual models for higher spin gauge fields [12, 13, 14 . 
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A Partitions and Discriminants 

A partition is a way of writing an integer n as a sumof positive integers where the order 
of the addends is not significant. The partitions of a number n correspond to the set of 
solutions (mi, ■ ■ ■ , m„) to the Diophantine equation 

1 mi + 2 1712 + ■ ■ ■ + n uin = n . (1 77) 

The number of solutions to eq. ( |177| ) is given by p{n) . This number is given by the generating 



function [M, 11 



i:Pin)x-=X{--— (178) 

n=0 n=l V-^ ) 

For the first values of n p{n) is given by 

p{n) = {1, 1, 2, 3, 5, 7, 11, 15, 22, 30, ■ ■ ■} . (179) 
For the first values of n the solutions to (|177|) are 



Ml = {{!)}, 

M2 = {(2,0), (0,1)}, 

Ms = {(3, 0, 0),(1, 1, 0),(0, 0, 1)}, 

M4 = {(4, 0, 0, 0), (2, 1, 0, 0), (0, 2, 0, 0), (1, 0, 1, 0), (0, 0, 0, 1)}, 
Ms = {(5, 0, 0, 0, 0), (3, 1, 0, 0, 0), (2, 0, 1, 0, 0), (1, 2, 0, 0, 0) 

(1, 0, 0, 1,0), (0, 1, 1,0, 0), (0, 0, 0, 0, 1)}, 
Me = {(6, 0, 0, 0, 0, 0), (4, 1, 0, 0, 0, 0), (3, 0, 1, 0, 0, 0), (2, 2, 0, 0, 0, 0), 
(2, 0, 0, 1, 0, 0), (1, 1, 1, 0, 0, 0), (1, 0, 0, 0, 1, 0), (0, 3, 0, 0, 0, 0), 
(0, 0, 2, 0, 0, 0), (0, 1, 0, 1, 0, 0), (0, 0, 0, 0, 0, 1)} . (180) 

Using this result we can now construct discriminants by 

p{n) n C_1^0■-l)■mfe,^ 
k=l j=l J ■""'kj- 

Permutations are obtained in a similar way 

Pin) n (_l\ij-iymkj 

^1 t^i j™'=^mfcj! 
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Grids for r = 4, n = 3 
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C Grids for r = 4, n = 4 
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Grids for r = 3, n = 4 







o- 


-o 






o- 


-o 















o- 


-o 


o- 

















f3 0\ 

3 

3 

VO 3/ 



/3 0\ 

3 

12 

VO 2 ly 

/2 1 ON^ 

12 

2 1 

VO 1 2/ 



16 



/3 0\ 

2 10 

2 1 

VO 1 2/ 





-o 


o- 


-o 




o- 


-o 





















o- 


-o 















o- 


-o 


o- 


-o 


o- 


-o 


o- 


-o 















1 









1 


1 


1 





24 





1 


1 


1 









1 


2/ 






1 









1 


1 


1 





48 








2 


1 




\o 


1 





2/ 


/I 


2 





o\ 




2 


1 
















1 


2 




VO 





2 















o- 


-o 


o- 


-o 











48 



/2 1 o^^ 

1110 

12 

VO 1 1 ly 

















-o 











^3 0>^ 

111 

111 

VO 1 1 ly 









-o 























o 






o 


-o 













o- 


-o 


o- 


-o 



















o- 


-o 




-o 








-o 























o- 













o- 


-o 


o- 


-o 










o- 




-o 


o- 




-o 











-12 



/3 0\ 

3 

2 1 

VO 1 2y 



24 



^3 0>^ 

2 10 

111 

VO 1 2y 



/2 


1 









1 


2 
















2 


1 




VO 





1 


2^ 






/^2 


1 





^\ 




1 


2 








12 








1 


2 




\o 





2 


1/ 




(3 








^\ 







2 


1 





48 








1 


2 




u 


1 


1 


1/ 






1 





o\ 







2 


1 





12 








2 


1 




u 








2y 


/I 


1 


1 


o\ 




1 


1 





1 




1 





1 


1 




VO 


1 


1 






/3 








ON 







1 


2 













1 


2 




VO 


2 










24 



/2 1 0\ 
12 

10 11 

Vo 1 2y 



56 



24 



24 



/2 1 ON^ 
1110 

111 

VO 1 2/ 

/2 1 ON^ 

111 

2 1 

VI 1 ly 





/I 


1 





1\ 







1 


2 





12 





1 


1 


1 




V2 








1/ 




/2 


1 





^\ 


48 





1 


1 


1 





1 


1 


1 




u 





1 


1/ 


/I 


1 





1^ 




1 


1 


1 










1 


1 


1 




Vi 





1 


l) 





48 



/2 1 0\ 

2 10 

12 

VI 1 ly 











o- 


-o 


o- 


-o 











-24 



-48 



/I 2 0\ 
1110 

12 

VI 1 ly 

/2 1 o\ 

12 

12 

VI 1 ly 



/I 2 0\ 

12 

12 

V2 ly 











o- 










(>• 















24 



/I 1 1 0\ 
1110 
12 

VI 1 ly 



57 



Grids for r = 6, n = 3 
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